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In numerous applications of energy harvesting via transformation of light into heat the focus recently 
shifted towards highly absorptive nanoplasmonic materials. It is currently established that noble 
metals-based absorptive plasmonic platforms deliver significant light-capturing capability and can be 
viewed as super-absorbers of optical radiation. Naturally, approaches to the direct experimental probing of 
macroscopic temperature increase resulting from these absorbers are welcomed. Here we derive a general 
quantitative method of characterizing heat- generating properties of optically absorptive layers via 
macroscopic thermal imaging. We further monitor macroscopic areas that are homogeneously heated by 
several degrees with nanostructures that occupy a mere 8% of the surface, leaving it essentially transparent 
and evidencing significant heat generation capability of nanoplasmon-enabled light capture. This has a 
direct bearing to a large number of applications where thermal management is crucial. 

Practical thermal platforms for solar energy harvesting and manipulation of thermal radiation are necessarily 
highly absorbing, but also affordable and scalable. Nanoplasmonic metamaterial absorbers 17 offer a route 
towards low material consumption and compact designs while maintaining large optical cross-sections 811 . 
They exploit the collective oscillations of electrons in metallic nanostructures - localized plasmons - for the highly 
efficient coupling of light. In the same time, individual plasmon nanostructures in these materials act collectively 
as an effective absorbing layer, i.e., a metasurface. The emerging field of thermoplasmonics then addresses 
applications of light-heated plasmonic nanostructures in photothermal therapeutics and drug release 11215 , ther- 
mal optical data storage 16 , enhanced catalysis 4,17,18 , magnetic recording 19,20 , solar thermal energy harvesting 6,2123 
and optoelectronics 24,25 . 

Whereas up to now the focus is on intricate design of nanoplasmonic metamaterials units, often by advanced 
electromagnetic simulations and highly demanding nanofabrication methods such as electron-beam lithography, 
recent advances in nanofabrication have led to simple and thereby cheap large-scale manufacturing methods of 
light absorbing plasmonic metasurfaces 26,27 . These approaches utilize amorphous arrangements of individual 
nanoplasmonic structures (meta-atoms) with simple geometries, which as a whole develop desired collective 
absorptive properties. When developing realistic light absorbing nanoplasmonic systems, optimization primarily 
focuses on the size, shape and material composition of the individual particles. While local plasmon-enabled 
thermal gradients, confined to the nanoscale, are well-documented in these systems 28 35 , large-scale (i. e., mac- 
roscopic) thermal effects due to plasmons are also of strong interest since their importance in e.g. catalysis has 
been demonstrated 36,37 . In the wake of the development of macroscopic plasmon metamaterial absorbers that aim, 
for example, at wafer-scale solar thermal harvesting, the comprehensive experimental characterization of such 
thermal effects comes handy. 

The prime aim when developing absorbing plasmonic architectures is to enhance the non-radiative decay of 
the excited nanoplasmons. This is typically done by reducing the size of plasmonic nanoparticles or by changing 
the dielectric function of the material 38 . Although noble metal nanoparticles of Au or Ag are widely used in 
plasmon metamaterials due to their well-pronounced resonances in the visible and near-IR, they might not be the 
best choice as absorbing platforms. Under white light or solar illumination plasmonic nanostructures having 
markedly lower but instead spectrally broader absorption may be equal to or even surpassing Au-based 
superabsorbers. 

Results 

We utilize simple large-area bottom-up nanofabrication method - hole mask colloidal lithography 27 - to prepare 
absorbing plasmonic metasurfaces, composed of amorphous arrays of Ni or Au nanodisks that are 20 nm high 
and nominally 110 nm in diameter (Figure la). Figure lb shows the absorption in these two types of nanoplas- 
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Figure 1 | (a) SEM micrograph of a plasmonic metasurface - amorphous arrangement of Ni nanodisks (scale bar - 5 um);(b) Absorption spectra of Au 
and Ni nanodisks and of carbon thin film of the same thickness as the height of the plasmon particles. Insets schematically show both types of 
absorbing platforms; (c) Experimental setup for measuring temperature and (d) a typical temperature response. 



monk metasurfaces, along with that of a carbon film of equal thick- 
ness that we use as the absorption benchmark. Au nanodisks layer 
shows strong and spectrally narrow absorption peak at the plasmon 
resonance. The Ni nanodisks system displays broader and generally 
lower absorption as a result of the shorter plasmon lifetime that is due 
to the strong electron-lattice coupling in Ni. 

Although optical spectroscopy gives straightforward quantifica- 
tion of light absorption in these metasurfaces, it understandably does 
not address the crucial information on the ability to generate elevated 
temperature at a macroscale due to such absorption. With the cur- 
rent study we put forward a new modality of characterizing mac- 
roscopic absorbing metasurfaces - that is, via their temperature 
increase upon illumination. We naturally consider that all absorbed 
light eventually transforms into heat in these systems. With this 
measuring modality we illuminate the surface with the solar simu- 
lator to benefit from potential omnichromatic absorption of the 
practical metamaterial absorber and track the resulting surface tem- 
perature with a thermal camera (Figure lc, further experimental 
details can be found in Methods section). Due to identical thermal 
radiation emissivity properties of the substrates (as opposed to the 
studied metasurfaces), we arrange the measuring system in a way that 
the ultra-thin absorbing plasmon metasurface is positioned at the 
backside of the sample and glass slide is facing light source and the 
thermal camera (c.f., Fig. lc). The temperature is further obtained by 
averaging the camera reading from the center part of the sample at 
which the illumination is homogeneous (typically, about 1 cm 2 ) first 
without the illumination (1 to 2 min), then with light on and reach- 
ing temperature saturation (about 3 min), and finally after light is off 
tracking the cooling-off kinetics (further 3 to 4 min). The result for 
the amorphous array of Au nanodisks as an absorber is shown in 
Figure Id, with vertical axis presenting the difference to the temper- 
ature in the dark (room temperature). At this point one clearly marks 



the macroscopic temperature effect due to the presence of plasmon 
metasurface and overall sample temperature increase up to 5-6 
degrees. We note here that temperature evolution in Fig. Id concerns 
the entire system - that is, plasmon metasurface, glass slide and the 
surrounding air altogether. Specifically, at the point of temperature 
saturation, a thermal equilibrium exists between these three media, 
and the temperature gradients are established within the glass sub- 
strate - the hottest being plasmon metasurface layer. 

To quantify the actual absorbed power we employ a simple numer- 
ical method. First we assume that, since it is thin and fully illumi- 
nated, the substrate has a homogeneous temperature. Any change of 
this temperature is directly proportional to a change of energy, per 
unit area, through the density of mass, specific heat, and thickness of 
the substrate. Differentiating this energy with respect to time gives 
the net thermal flux, heating (or cooling) the substrate. In principle 
the substrate has a constant positive (i.e. heating up) flux of heat 
coming from the absorbing metasurface upon solar illumination. A 
negative heat flux is present when the substrate has a higher temper- 
ature than its surroundings. This negative flux is proportional to the 
temperature difference between the sample and its surroundings 
through two main terms: convective and radiative heat transfers. 
The convective term is a function of the temperatures to the power 
of unity, the radiative - to the power of four. However, since we are 
interested only in a small temperature window far from absolute 
zero, the radiative term can also be assumed to be linear with respect 
to the temperature inside this measurement window. Now, plotting 
the thermal flux against the substrate temperature displays two 
regions of interest (Figure 2a, left): the data corresponding to heating 
and cooling of the substrate, respectively. Calculating the offset 
between linear fits to the heating and cooling data earns the heating 
power. Figure 2a schematically outlines the general procedure of 
extracting the absorbed power with the given system (the actual fit 
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Figure 2 | (a) Schematic outline of the absorbed power measurement with two complimentary methods - temperature reading, fitted with numerical 
model (left) and optical absorption, combined with the emission spectrum of the light source (right). Net heat flux in (a) corresponds to the time- 
differentiated temperature in the blue and red highlighted rectangular regions of the temperature-vs.-time plot, multiplied by the factor for temperature- 
to-energy conversion, (b) Absorbed power from thermal (light grey) and optical (dark grey) measurements (error bars shown for thermal 
measurements), green bars mark the ratio between the two (right axis). 



to experimental data is shown, Figure 2a, left) and cross-referencing 
it to the measured optical absorption (Figure 2a, right). The absorbed 
powers for three systems, extracted both from the temperature data 
(light grey) and optical absorption (dark grey), are presented in 
Figure 2b. Here we compare two nanoplasmonic platforms (amorph- 
ous arrays of Au or Ni nanodisks) and carbon film. For completeness 
we also plot the ratio of the two (green bars, with the corresponding 
vertical axis to the right). The latter visually confirms the consistency 
between the thermally derived absorption and the optically measured 
one in all three systems. In order to keep the generality of the 
approach and maintaining the robustness of the absorbed power 
extraction against the particularities of the exact solar thermal system 
under consideration (i.e., exact substrate thickness, absorbing meta- 
materials composition etc.), we rely on the best fit of the experimental 
temperature data that the model provides. This, in general, under- 
estimates the extracted absorbed power as verified by the direct 
optical absorption measurements (c.f Fig. 2b). Furthermore, the ratio 
of thermal-to-optical retrieved absorption power varies between the 
samples. This variation is largely due to uncertainties in the optical 
absorption measurements. The optical absorption is measured with 
an illumination spot size of only a few mm 2 . Slight variations in 
particle coverage over the sample surfaces therefore have significant 
impact on the measurement result. In contrast, the thermal method 
averages absorption over the larger area, which makes it less sensitive 
to such inhomogeneity. 



We further establish that the quantitative aspects of the absorptive 
properties of a given thin functional layer, in principle, of unknown 
composition and morphology can be directly extracted from the 
thermal data. For example, comparing absorptive Au metamaterial 
surface with carbon film, one concludes that the carbon film absorbs 
4.5 times more than Au metamaterial. In our case this is also directly 
verified by the measured optical absorbance for all three systems. 
Excitingly, if the morphology of the absorptive layer is known (i.e., 
film thickness and/or nanostructures size and coverage, so the overall 
content of the absorptive material is quantified) - straightforward 
comparison of the heat generation ability of these ultrathin mac- 
roscopic functional layers can be made. In our case, normalizing 
the thermal power on materials coverage both for nanoplasmonic 
metasurfaces and carbon film, we estimate nanoplasmonic systems 
to deliver 3.4 times higher absorption/thermal power per surface unit 
given that they have a surface coverage of 8 per cent. 

Now that we outlined the method to correlate total broadband 
absorption and thermal effects, we turn back to the plasmon-enabled 
heat-generation platforms to mark the possible route towards an 
optimal plasmonic 'solar heater'. As mentioned earlier, gold- and 
silver-based plasmonic nanoarchitectures typically display strong 
and spectrally well-defined region of light absorption that is linked 
to the plasmon resonance. By contrast, metals like Ni present mark- 
edly broader plasmon feature in absorption due to higher ohmic 
losses. We use the latter property to design a simple plasmonic solar 



SCIENTIFIC REPORTS | 4 : 51 1 1 | DOI: 1 0. 1 038/srep05 1 1 1 



3 



temperature system that favorably compares to structurally identical 
one, but made of Au. Namely, we use presented method to compare 
the absorbed power in amorphous arrays of Au and Ni nanoellipses 
of three sizes: (i) 80 X 60 nm, (ii) 112 X 60 nm and (iii) 140 X 
100 nm. Figures 3a-c compare longitudinal and transversal optical 
absorption of Au (right panels) and Ni (left panels) nanoellipses. (the 
overall array morphology for both materials is depicted in Fig. 3d). 
Even though both longitudinal and transversal peak optical absorp- 
tion in Ni systems are visibly lower (compare left and right panels for 
all nanoellipses sizes), significant spectral broadening of the absorp- 
tion in Ni makes total areas under the peaks, which in this sense 
correspond to the actual absorption over the entire solar spectral 
range, noticeably larger (c.f. the solar irradiance schematics in 
Figs. 3 a-c at the background). In general, larger plasmonic particles 
are known to have larger scattering cross-sections than smaller 
ones 39 , thereby loosing in their absorption cross-section. Due to their 
strongly damped plasmon resonances, Ni nanostructures deviate 
from this generic trend. This fact directly reports in the thermally 
measured data (Figure 3e), where upon non-polarized solar simu- 
lator illumination 100 X 140 nm Ni nanoellipses generate 30% more 
heat than corresponding nanostructures of Au. Note that thermal- 
to-optical absorbed power ratio tends to be the same as in the case of 
Ni and Au nanodisks and the carbon film (compare the thermal-to- 
optical ratios in green on Figs. 2b and 3e). Overall, the comparison of 
differently shaped and sized nanostructures of various materials sig- 
nals that Ni-based plasmonic absorbers would potentially deliver 
appreciably large and spectrally broad heat- generation capability in 
practical applications. 

Efficiency under oblique incident light is another practically 
important property of a solar absorber. In Figure 4 we present gen- 
erated heat from nickel nanodisks and nanoellipses, illuminated with 
light incident at 0, 30 and 60 degrees off normal, with the plasmon 
metasurface at the backside (i.e., in the same measuring geometry as 
in all discussed cases here, see Figure lc). Figure 4 presents data for 
nanoellipses oriented both with long and short axes in the plane of 
incidence. Interestingly, the efficiencies of the thermal power gen- 
eration by the plasmon metasurfaces at oblique illumination are 
slightly higher than what can be expected from the regular photo- 
voltaics cells in the similar light conditions. In fact, in photovoltaics 
the decrease in insolation, with a cosine dependence on the angle of 
incidence, is one of the main causes for the lowered resulting effi- 
ciencies 40 . Cosine dependence dictates an insolation loss of 50 per 
cent at 60 degrees angle of incidence from normal. Plasmon meta- 
surfaces (both of nanodisks and nanoellipses) display about 30% 
decrease in the absorbed power generation when the illuminating 
light changes incidence from 0 to 60 degrees from normal, thus 
positively comparing with the relevant numbers in photovoltaics. 
It would seem that the increased optical path upon the grazing incid- 
ence, possibly combined with the waveguiding modes in the glass 
substrate, positively contributes to the absorptive ability of plasmon 
thermal metasurfaces. 

These results bring forward the potential of other 'non-conven- 
tional' plasmonic materials - such as Pd, Pt, Co and others - as 
plasmon-enabled thermal metamaterial platforms for solar energy 
harvesting 41,42 . They typically feature very moderate absolute absorp- 
tion in a broad spectral range (c.f., Fig. lb) that also comfortably 
covers solar radiation spectrum. This leaves macroscopic areas that 
are functionalized with thermal plasmon nanostructures essentially 
transparent. In the same time, these nanostructures provide superior 
heat-generation ability. With the simultaneous access to design 
and affordable fabrication of truly large-scale thermal plasmonic 
metasurfaces from these materials, the concepts of plasmon-enabled 
heating might spill into applications of thermal surfaces, for example, 
in energy-efficient technologies for construction and industrial 
design. Another exciting prospect is in thermal plasmonics of 
nanoferromagnets for the manipulation of magnetic order. Indeed, 



plasmon-induced temperature effects are currently actively explored 
with magnetoplasmonics 1 '' 20 With the current work we show 
that Ni as efficient nanoplasmon heat antenna 43,44 has substantial 
heat-generation ability and thus might find direct use in promoting 
strong nanoscale temperature gradients for effective magnetization 
manipulation. 

Discussion 

Summarizing, we devised a modality to comprehensively character- 
ize the absorbing and heat-generating ability of various ultra-thin 
absorbers, most importantly including plasmon-enabled mac- 
roscopic absorbing metasurfaces. It is also very appealing that the 
method allows for quantitative comparison of fundamentally differ- 
ent absorbing systems. Further, we used the method to mark the 
possible route towards potentially high-efficient omnichromatic 
and omnidirectional plasmon absorbers that utilize Ni rather than 
Au. It is also appealing to exploit high spectral tunability of various 
plasmonic nanostructures to address various spectral regions for the 
most efficient heat power generation with such metasurfaces. There 
is a wide and immediate implication of these findings - that is, this 
work provides very simple tool to increase the credibility of the 
studies addressing super-absorbing platforms - metallic, dielectric 
or semiconducting. It aids to make a direct step from claiming 
highly absorptive properties to the actual proof of the resulting 
increased temperature. Specifically, studies using a spectrally variable 
monochromatic light source instead of a solar simulator would in 
addition allow wavelength-resolved measurements of the actual heat 
generation. 

All these powerful features open up new possibilities for the reli- 
able evaluation and practical implementation of ultra-thin metama- 
terial superabsorbers in photocatalysis, energy- efficient design, 
thermal magnetization control and other applications where thermal 
management is important. 

Methods 

Plasmonic nanodisks and nanoellipses were prepared with the hole-mask colloidal 
lithography process as it is described elsewhere 27 . Polysterene particles with 110 nm 
diameter were used to fabricate the masks for the disks. 60, 80 and 100 nm particles 
were used for the masks of 84 X 60, 112 X 80 and 140 X 100 nm ellipses respectively. 
The thickness of all structures was 20 nm. 0.21 mm thick circular borosilicate cover 
glasses with 25 mm diameter (D 263 M, Menzel-Glaser) were used as substrates for 
the disks and 0.4 mm thick 15X15 mm cover glasses of the same type were used for 
the ellipses. The samples were entirely homogeneously covered by plasmonic nano- 
particles. Samples with carbon thin films were prepared by e-beam evaporation of 
carbon on 0.5 mm thick fused silica substrates followed by annealing in Ar atmo- 
sphere at 800°C for 10 min. 

Optical absorption measurements were done using Cary 5000 spectrophotometer 
[Varian] with an integrating sphere [DRA 2500, Varian] equipped with the small- 
spot kit. The samples were either positioned using the center mount sample holder or 
placed at the reflectance- and transmission ports alternately. In both cases specularly 
reflected light was included. The baselines were recorded with an empty sample 
holder in place. SEM micrographs of the plasmonic nanostructures were obtained 
from the same samples as the other measurements were made on using a Quanta 200 
FEG ESEM, FEI. Particle sizes and number densities were determined from the 
micrographs using the software Image J 1.44o. 

Time resolved photo- induced surface temperatures of the samples were measured 
with a setup consisting of a solar simulator [ssl50, Sciencetech - see the output 
spectrum in Appendix] as light source, thermal cameras A20, FLIR. (for nanodisks 
samples and carbon films) and A645 sc, FLIR (used for nanoellipse) were used to 
monitor the sample temperature and a sample holder with the sample. The samples 
were positioned in the center of the solar simulator light beam at normal angle to 
ensure homogeneous illumination in the center part of the samples. The intensity of 
light was calibrated to 2 suns (2000 W/m 2 ) using a photodiode [Sciencetech] . The 
absorbing metamaterial side of the sample (i.e. the backside) was faced away from the 
light source, and the thermal camera was directed towards the front side of the sample 
at an angle large enough for the IR-radiation specularly reflected by the sample into 
the camera to be that of the ambient environment and not from the camera itself or 
the solar simulator. The ambient environment has a known and constant emission as 
opposed to the thermal camera or the solar simulator, thus allowing reliable thermal 
measurements. Measuring the temperature of the samples at an angle does not affect 
the reading since both the thermal emission and projected area have a cosine 
dependence on the angle of observation. The thermal camera in principle calculates 
the temperature from measured radiosity (radiated power per unit area) and since 
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Figure 3 | Longitudinal (blue) and transversal (red) optical absorption in Ni (left) and Au (right) ellipses of sizes (a, i) 60 X 84 nm, (b, ii) 80 X 112 nm 
and (c, iii) 100 X 140 nm. nanoellipses. Schematics of the mutual orientation of nanoellipses axes and the polarization of the incoming light in absorption 
measurements. Solar irradiance (not to scale) is shown schematically in the background as a guide, (d) SEM image of a typical nanoellipses arrangement, 
prepared by HCL. (e) Thermally measured absorbed power in Au (orange) andNi (grey). Ratios to spectroscopically-retrieved absorbed power are shown 
with green bars (right axis). Error bars on absorption data indicate the standard deviation of at least 3 measurements on the same sample. Error 
bars on thermal-to-optical absorbed power ratios indicate 9% deviation from the weighted average. 
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— nanodisks 

— nanoellipses, long axis in incidence plane 

— nanoellipses, short axis in incidence plane 

Figure 4 | Absorbed power (measured with the thermal camera) for Ni 
nanodisks (blue) and nanoellipses (red and green) with light incident at 0, 
30 and 60 degrees from normal. Nanoellipses are oriented with long 
(green) or short (red) axes in the plane of incidence. 

both the perceived area and the perceived radiosity scales equally with respect to 
observing angle the both effects cancel each other. 

The emissivity of the sample front side (glass surface) was verified by heating the 
substrate above 50 J C and measuring its temperature using a thin k-type ther- 
mocouple with its hot junction adhered to the substrate by Ag-glue and its cold 
junction kept at a controlled 0°C. At the same time the temperature of the substrate 
was measured with the thermal camera. As the temperature readings were equal, the 
correct emissivity can be recorded in the real experiments. To reduce the heat dis- 
sipated from the sample via thermal contacts, the sample holder balanced the sample 
between two razorblade edges and a sharpened wire tip. 

The thermal power delivered to the substrate was estimated using a MatLab script. 
Curve fitting was done using the MatLab internal function 'fit*. Material properties 
consisted of the density of mass and the specific heat capacity. Values for the cover 
glasses were provided by the manufacturer. Fused silica values were taken from the 
literature 45 . 
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